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ABSTRACT 

A Field Programmable Gate Array implementation of the Gabor-type filter is presented. The implementation 
uses the forward Euler approximation. In, this paper the details of the realization of Gabor filter architecture 
using the floating - point operations are discussed The paper presents, a brief mathematical overview of the 
Gabor-type filters which involves floating point arithmetic operation analysis results, the architectural details 
of the Gabor-type filter, and, the structure of the FPGA implementation is given. Finally implementation 
results are discussed. Simulation and Synthesis is done using Xilinx ISE design suite. Verilog HDL will be 
used as a description language for mapping algorithm in VLSI and hardware implementation on SPARTAN -3 E 
FPGA. 
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I. INTRODUCTION 

In image processing, a Gabor filter, named after Demis Gabor, is a linear filter used for the edge 
detection. Frequency and orientation representation of Gabor filters are similar to those of the human visual 
system, and they have been found to be particularly appropriate for texture representation and discrimination. 
In the literature, FIR [3] -[6] and CNN - based [7] digital implementation methods for Gabor filters are 
available, where [3]-[6] demonstrated graphics processing unit (GPU) and application - specific integrated 
circuit implementations respectively However there are some difficulties in the realization of FIR Gabor-filters 
using fixed-point operations, such as accuracy, large template requirements for high quality factor etc. The 
difficulties have been overcome by Gabor-type filters realized with CNN using floating point operations, which 
has similar properties to Gabor-filters. This paper deals with mathematical overview of the Gabor-type filter, 
which involves floating point arithmetic operation analysis, results and the architectural details of the Gabor- 
type filter are presented and the structure of the FPGA implementation is given. Finally implementation results 
are discussed 

II. CNN - BASED GTF ARCHITECTURE 

A DT-CNN-based GTF structure was proposed in [2] as an implementation of [1] whose signal-flow 
diagram is shown in figure 1 . below. 

Where, lcx,. D arid co^- are die center frequencies ofx and y axes respectively 
JR= v^vel ength of die sinusoidal factor. 
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Fig. 1 DT-CNN-based GTF structure signal flow graph [22] 



THE IJES 



www.theijes.com 
Page 25 



The IJES 



Design And Implementation Of Gabor... 



The signal flow graph of the above figure 1 is as shown as a logical block diagram of one CNN - based GTF 
computation unit in figure 2. A DT-CNN -based GTF structure forms the basis for the logical block diagram of 
figure 2 which is used in this paper for the design and implementation of the GTF architecture. 
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Fig. 2. Logical block diagram of one CNN-based GTF computation unit 

Two different approaches are used in the logical realization of the fig 2: direct implementation and resource 
sharing. In the former, a computation tree is implemented directly and, in the latter, considering the similarity 
of the two diagrams; one computation tree is realized and multiplexed for the computation of real and 
imaginary parts (fig. 2) 

The multiplexing halves the number of multiplier and adder resources; however, twice the pixel frequency is 
required for computation. The input bu^ is also multiplexed with zero, as it is only required for the calculation 
of the real part. All adders and multipliers are registered to form a pipeline. The latency of the processor is 
three pixel clock cycles as there are six pipe stages. Finally, the real and the imaginary parts of the resulting 
values are demultiplexed to their respective output registers. 

III. RESULTS AND DISSCUSSION 
3.1 SIMULATION RESULTS 

The architecture was written with Verilog HDL and synthesized by Xilinx ISE 14.2. After 
synthesizing this Verilog HDL code, tests were done with MODELSIM. The following are the simulation 
results obtained. 

3.1.1 Simulation and synthesis results for select line's' zero 

The below window represents the input values for select line zero. The calculated floating-point 
binary numbers for each multiplexer output is shown. Thus indicating, for the GTF architecture the results are 
accurate, matching both the manual and the synthesis results exactly. 

Fig3. by forcing the clock as ' 1 ', period of ' 100|is' and select line 's' zero, the results of the multiplexer and 
de-multiplexer pipeline stages of the GTF architecture are obtained representing the real part of the 

architectural design. 
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Fig. 3 representing input values for the real part of the system with select line zero 
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Fig 4 represents the boundary conditions alpha, beta, lambda values of real term with select line's' zero. 
Also resulted multiplexed registered outputs and the demultiplexed real term [15:0] is obtained 
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Fig. 4 representing the boundary conditions and registered outputs [15:0] real term of the system with select 

line zero 

3.1.2 Synthesis and simulation results for select line's' one: This represents the imaginary value of the 
demultiplexed output. Below window represents the input values for the imaginary term of the multiplexer and 
de-multiplexer. 

Fig. 5 : by forcing the clock as '1', period of ' 100|as' and select line 's' one, the results of the multiplexer 
and de-multiplexer pipeline stages of the proposed architecture are obtained representing the imaginary part of 

the architectural design. 
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Fig 5 representing the real part of the system with select line one 



Fig 6 represents the boundary conditions alpha, beta, lambda values of imaginary term with select line's' 
one. Also resulted multiplexed registered outputs and the demultiplexed real term [15:0] are obtained. 
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Fig. 6 represents the boundary conditions alpha, beta, lambda and imaginary term 



3.2. FPGA IMPLEMENTATION RESULTS 
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In order to form Gabor type filter structure shown in fig 2, multipliers and adders are needed. The 
structure is formed by blocks shown in fig 1. The logical block diagram of fig 1 is depicted in fig 2 This 
architecture was written with VHDL and synthesized by Xilinx ISE 14.2. After synthesizing this VHDL code, 
tests were done with MODELSIM. When this process unit is synthesized, placed and routed for XC3S400 
FPGA, a 50-MHz clock rate can be achieved. The architecture can cover 186 (1%) slice flip flops and four 
mult 18x18 (20%) on FPGA. 




Fig. 7 FPGA implementation of GTF architecture with the hardware resulting in real part of the demultiplexer 

with select line's' zero. 




Fig. 8 FPGA implementation of GTF architecture with the hardware resulting in imaginary part of the 

demultiplexer with select line 's'one. 
The fig. 9 below shows the top module of the system architecture. 
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led test 




led test 



Fig. 9 Top module technology schematic 
The Fig. 10 below shows the internal architecture of the top module 




Fig. 10 The internal architecture of top module 
The fig. 1 1 shows the schematic diagram of the GTF architecture. 
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Fig. 1 1 The schematic diagram of the GFT architecture 



The fig. 12 shows the power analysis results for the desugned architecture. 
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Fig 12 Power analysis 



IV. CONCLUSION 

In this paper, a brief mathematical overview of the Gabor-type filters is presented; floating point 
arithmetic operation analysis results and the architectural details of the Gabor-type filter are given. The 
structure of the FPGA implementation is given. Finally implementation results are discussed. A CNN- Gabor 
type filter realization method was proposed coded in VHDL and the simulation results are obtained resulting in 
accuracy. 

Digital implementation of Gabor-type filters was proposed, which drastically reduces the number of multipliers 
required for FPGA implementation. The Gabor filters are suitable for a real-time application such as optical 
character recognition (OCR), facial recognition or license plate recognition and any other system that requires 
two-dimensional band-pass filters. 
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